Spin Glass Ordering in Diluted Magnetic Semiconductors: a Monte Carlo Study 
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We study the temperature-dilution phase diagram of a site-diluted Heisenberg anti-ferromagnet 
on a fee lattice, with and without the Dzyaloshinskii-Moriya anisotropic term, fixed to realis- 
tic microscopic parameters for IIBi-^ Mn^Te {IIB=Cd, Hg, Zn). We show that the dipolar 
Dzyaloshinskii-Moriya anisotropy induces a finite-temperature phase transition to a spin glass phase, 
at dilutions larger than 80%. The resulting probability distribution of the order parameter P{q) is 
similar to the one found in the cubic lattice Edwards- Anderson Ising model. The critical exponents 
undergo large finite size corrections, but tend to values similar to the ones of the Edwards- Anderson- 
Ising model. 

PACS numbers;75.10.Nr, 75.40. Mg, 75.40. Gb, 64.60. Cn, 64.60.Fr 



Although most theoretical investigations on spin glass 
(SG) systems have focused on models based on 
Ising spins, systems that have been investigated in 
the recent period, like the diluted magnetic semi- 
conductors IIBMnTe series, JIB — Cd,Hg,Zn, 
i.e. Cdi.^Mn^Te |,|, Hg^^^Mn^Te § and 
Zrii-x Mrix Te m (and somehow even typical experi- 
mental samples pi like CuMn, AgMn, Eu^ Sri-^ S), 
are closer in nature to continuous Heisenberg spins. Early 
computer simulations suggested that in three spatial di- 
mensions neither systems with local interactions and 
Heisenberg or XY |^,^ spins, nor systems with long- 
range RKKY interactions |lO| undergo a finite temper- 
ature SG phase transition. This fact could be poten- 
tially annoying from a phcnomenological point of view, 
but it becomes acceptable after noticing that a small 
anisotropic interaction, neglected in the above calcula- 
tions, could induce a finite temperature SG phase tran- 
sition. 

Still the situation is not crystal clear: the work of ||Tl| ] 
claimed that the most important anisotropic coupling 
in II B Mn Te materials , the Dyalozhinskii-Moriya 
(DM) interaction, is not able to induce a SG phase. On 
the contrary a theoretical analysis of experimental 
data on the II B Mn Te series was used to suggest the 
presence of SG ordering in 31? Heisenberg spin glasses 
(for finite temperature and no anisotropics): recent nu- 
merical simulations |Q| support the existence of a chi- 
ral phase transition in such systems. The role of the 
anisotropy was reconsidered in [T5| ], where the Heisen- 
berg spin Edwards- Anderson (HEA) model was consid- 
ered with the addition of a random pseudo-dipolar inter- 
action: clear signatures of a finite temperature SG phase 
were found. This result is not consistent with the one 
of ijllj (that was considered as being based on a realistic 
modelization of IIB MnTe, even if the direction of the 
DM vectors, see was chosen at random, while they 
should be periodic along the lattice Q). We remind at 



last that a diluted Ising anti-ferromagnet on a fee lat- 
tice has been studied in ||l6|, where the signature of a 
SG phase transition for low enough densities has been 
detected. 

We have taken here the point of view of trying to be as 
realistic as possible, analyzing a model as close as possible 
to the experimental samples. We show that non-random 
DM terms (selecting a realistic value for the anisotropy) 
are able to induce a SG phase transition in a 3D Heisen- 
berg spin glass on a fee lattice. We analyze and discuss 
in detail the values of critical exponents: the experimen- 
tal results for the IIBi-x Mux Te materials [||^] are in 
good agreement with the most accurate calculations for 
the Edwards- Anderson model with Ising spins (lEA) on 
the cubic lattice (j^ = 1.8 ± 0_^, ?? = -0.26 ± 0.04), 
but the numerical simulations of [|l6| and of yielded 
ly « 1.0. We will show that the numerical calculation of 
the critical exponents on the accessible lattice sizes suffer 
from serious finite-size corrections, and that a systematic 
analysis of the numerical data establishes clear trend to- 
wards values of larger than the ones found in |]l5| , |l6| , 
and close to the experimental values. 

The site-diluted anti-ferromagnetic (AFM) Heisen- 
berg model on the fee lattice, with and without DM 
anisotropy, is a model for the IIBi-x Mux Te series, 
where the Mn atoms form an fee lattice with local- 
ized (Heisenberg) spins interacting through short-range 
(super-exchange) AFM terms, while the magnetically in- 
ert IIB atoms randomly replaces the Mn over the lattice. 
An AFM interaction on the fee lattice is frustrated, and 
gives rise to some interesting order-disorder phenomena 
||l8|, both with Heisenberg Q and Ising spins The 
dilution disorder deletes some of the sites on the system, 
thus providing the random combination of frustrated and 
unfrustrated plaquettes, that is believed to be essential 
for SG ordering. The Hamiltonian of the system is 
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where the fields S — {Si, S2, S3), Si, S2 and S3 real 
with = 1, represent the spin of the Mn atoms, and 
J > 0. A lattice site is randomly occupied by a spin 
with probability p. The sum labeled by {x, y)' runs over 
the pairs of occupied nearest neighboring sites of the 
lattice. The unit-length vectors Rx-y specify the DM 



anisotropy, and they verify Rx % 
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||12[ we set R,^ ^ n) = i~7^^ — 7?'0)' while the other 
five independent vectors are obtained using the three-fold 
rotation symmetries of the lattice. In Ref. the ratio 
D/J has been estimated to be 0.054 for Zn.77Mn.33re 
and Cc?. 77 Mn. 33 Te, and is very mildly dependent on the 
composition of the sample: we have fixed it to 0.06 for 
simplicity. As the local magnetic field acting on spin Sx 
is (see (0)) 
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it is easy to implement a heat-bath algorithm and the 
over-relaxed micro-canonical algorithm of [|l9|. We have 
found that the combination of these two updates tremen- 
dously reduces the thermalization effort. In the pro- 
duction runs we have performed a full-lattice heat-bath 
sweep followed by 19 over-relaxed updates, that will be 
referred in the following as an elementary Monte-Carlo 
step (EMCS). In order to define the observables, it is 
useful to consider a replica (i.e. a thermally independent 
system with the same set of occupied sites, that we de- 
note S{x)). The measured observables can be most eas- 
ily described in terms of the following three basic fields: 
the tensor field (TaAx) = S'^S^ - iJ"-'', if the lat- 
tice site X is occupied, and zero otherwise), the tensorial 
overlap {Oa.pix) = S^S^), and the scalar overlap field 

iq{x) = Sx-~Sx). 

The rationale for studying the tensor field t{x) j2^ is 
that previous studies of AFM diluted systems on the fee 
lattice |l^], showed that only for moderate dilutions the 
system ceases to develop AFM ordering. The tensorial 
magnetization is an ideal order parameter to check this 
possibility, since it will be non- vanishing for any conceiv- 
able type of AFM or helicoidal ordering. Moreover, it 
would also work on more sophisticated, yet trivial situ- 
ations like the found in D = 0, p = 1 jl^. The tenso- 
rial overlap p2| is most adequate to study the isotropic 
(D — 0) case, since in this situation the Hamiltonian 
posses a 0(3) global symmetry: when the anisotropy is 
switched on the symmetry reduces to Z2, and the use of 
the scalar overlap becomes natural. For all three fields, 
one can define straightforwardly the corresponding sus- 
ceptibility. Binder parameter and a finite lattice correla- 
tion length raJlTpl]. 
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FIG. 1. Tensorial {Bt) and scalar (Bg) overlap Binder pa- 
rameter versus temperature for sizes L = 8, 16 and dilution 
p = 0.9 (upper part) and p = 0.7 (lower part). 

The model (|l|) without impurities {p = 1.0, D = 0.06 J) 
undergoes a phase transition at Tc{p = 1) « 0.60J from 
a paramagnetic phase to an AFM phase, as shown by the 
behavior of the tensorial magnetization. For larger dilu- 
tions a lower temperature value needs to be reached in 
order to exit from the paramagnetic phase: the critical 
line, Tc{p), will eventually reach zero temperature at the 
percolation threshold for the magnetic ions {pc ~ 0.2). 
The first question is for which dilution the system forgets 
its global AFM ordering. In order to answer this question 
we have performed slow annealings in 60 samples (and its 
corresponding rephcas) at dilutions p — 1.0,0.9 and 0.8, 
in lattices L = 8 and 16; at p = 0.7 and p = 0.6 (that will 
be shown to be in the SG compositional range), we have 
annealed 700 samples. The results for the Binder cumu- 
lant of the tensor and scalar overlap fields at p = 0.7 and 
p = 0.9 are displayed in figure ^ At p = 0.9 for both ob- 
servables we find a low temperature, AFM ordered phase, 
since the tensorial magnetization is non- vanishing. There 
is a strong dip close to the phase transition point, which 
probably is very plausibly of first order. On the con- 
trary for p — 0.7 it is clear that the tensorial magneti- 
zation is no longer an appropriate order parameter. For 
p = 0.8 (not shown in the plot), our results indicate a 
cross-over regime between the two situations. Therefore, 
the low temperature phase turns from AFM to SG at 
0.7 < Pc < 0.8, similarly to what happens in the Ising 
case |l|]. Also for L>^0,p< 0.8 we do not find an AFM 
ordered phase. 

In order to quantify how strong the effect of the 
anisotropy is, we compare the system at p = 0.7 with 
D — and D — 0.06 J. In figure || we show the correla- 
tion length of the tensorial overlap in units of the lattice 
size. This operator should be zero in the paramagnetic 
phase, diverging in the SG phase, and at the critical point 
reaches a finite universal value. In the isotropic case (see 
zoom in upper part of figure |^) the crossing point of 
the L = 12 and 16 lattices is not clearly resolved, and 



2 



their respective crossings with the L = 8 curve shifts to 
lower temperature with growing lattice size. Moreover, 
the data (not shown in the plot) for the Binder Cumu- 
lant of the tensorial overlap rapidly grow at the crossing 
temperature of the correlation length, but then saturate 
at a value which decreases with the lattice size, without 
a crossing, similarly to the results shown in [ p^ , where it 
has been shown that the chiral-glass phase appears pre- 
cisely at the temperatures at which the Binder cumulant 
of the tensorial overlap grows. 




FIG. 2. Correlation length of tensorial overlap (in units of 
lattice size) vs. temperature for D — (upper part and zoom) 
and in presence of anisotropy D = 0.06 J (lower part). 

On the other hand, with a 6% DM anisotropy (see 
the lower frame of figure ||), we find a neat crossing of 
the correlation length (the tensorial overlap Binder cu- 
mulant has a marked dip, in contrast with the scalar 
overlap shown in the lower frame of figure 1). As the 
phase transition for the anisotropic system occurs at a 
temperature 80% higher than the one close to the cross- 
ings of the D = case, that according to Q signal a 
real chiral-glass phase transition, the natural conclusion 
is that the DM anisotropy is not a smooth perturbation 
that reveals a hidden chiral-glass ordering. 




FIG. 3. Probability distribution function of scalar over- 
lap for L = 8, 12, 16 at p = 0.6, T = J/4.5 ^ 0.78Tc and 
D = 0.06J. 



In order to characterize more precisely the SG phase 
we have studied the distribution of the scalar overlap at 
p = 0.6, D = 0.06 J and T = J/4.5 w O.TST^. At this 
temperature we have estimated the mean thermalization 
time in the L = 16 lattice, by considering a logarithmic 
plot of the mean overlap susceptibility of 64 samples, as 
a function of MC-time, starting from a random configu- 
ration, and we have found it to be of order 250 EMCS. 
After that we have performed a run with 800 samples, 
with L = 8, 12, 16, performing 8000 EMCS on each sam- 
ple, and taking a measure every 4 EMCS. We display 
P{q) in figure ^ is central part is remarkably stable for 
growing lattice size. Therefore, on the lattice sizes that 
we are able to thermalize, the pattern we obtain is com- 
pletely analogous to the one found for the Ising EA model 
in 3D m. 

Due to the global Z2 symmetry of the Hamiltonian 
(|l|), one would expect it to belong to the same univer- 
sality class of the lEA model in 3D, which seems even 
more plausible from our measures of the P{q)- To further 
investigate this relation we have measured the critical ex- 
ponents, in the dilution range where we definitively find 
SG ordering, namely p = 0.7 and 0.6. Since we have 
at our disposal only a narrow range of lattice sizes, it is 
important to use a finite size scaling analysis that allows 
to study the scaling corrections. We have used the quo- 
tient method of [|l|, that has been particularly useful 
in the study of scaling corrections in disordered systems 
p^ . We measure an operator O, diverging in the infi- 
nite size limit at criticality as |T — Tc|~^°, on two finite 
lattices of sides L and sL, and we select the temperature 
value where the two correlation lengths in units of the 
lattice size coincide (see the crossing of figure ||) . For the 
quotient of these two measures we have 
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where a; > is related to the first irrelevant operator 
in the Renormalization Group sense. The main advan- 
tage of the relation (||) is that the large statistical cor- 
relation between the measurements of O and ^ allows 
to measure the quotient with sufhcient accuracy as to 
uncover the scaling corrections. In our Z2 symmetric 
case we have of course used the scalar overlap correla- 
tion length. Our results are displayed in table |, and 
they do show the presence of significant scaling correc- 
tions {v is computed using the temperature derivative 
of ^g, 77 from the susceptibility of the scalar overlap). 
Since we only have few lattice sizes it is meaningless to 
try an infinite size extrapolation as the one of |^^. We 
should still mention that the critical exponents obtained 
by a simple log-log fit (for instance, with data measured 
at the maximum of the specific-heat) is roughly equiva- 
lent to an average of the transient exponents displayed 
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in our table. Therefore the value of « 1 found with 
Ising spins | Iqj or in the HEA model with pseudo-dipolar 
anisotropy la, is most probably a preasymptotic value. 
In fact, the best available results for the IE A model in 
3L» |l|], V = 1.8(2), ?7 = -0.26(4) are plausible infinite 
volume extrapolations for our results. However, in order 
to definitively elucidate this point, it would be helpful 
to use the extrapolation method of [^,0, that allows 
to work in the paramagnetic region with a significantly 
smaller thermalization effort. 

We have shown for the first time that dipolar DM 
anisotropic local interactions are able to induce a SG 
phase transition for Heisenberg spins in three dimensions. 
This result has been obtained with the very small re- 
alistic value of the anisotropy coupling constant in the 
IIB MnTe series. Given the dramatic effect of this 
small perturbation term in the Hamiltonian, we suggest 
that the chiral-glass mechanism proposed in [|l3| is over- 
whelmed by the neglected DM term. We have studied the 
temperature-dilution phase diagram of the Hamiltonian 
(|l|) in a large dilution range. We have found that the 
low-temperature phase changes from AFM to SG order 
between p = 0.8 and p = 0.7. We have used a com- 
bination of micro-canonical and heat-bath Monte Carlo 
update, that have allowed us to thermalize a L = 16 lat- 
tice at T = 0.78Tc, in the SG phase. We have measured 
the distribution of the overlap, finding results analogous 
to the ones of Ising EA 3-D model on similar lattice sizes. 
We have given an estimate the critical-exponents on the 
SG dilution range. Our results suffer from severe finite 
size corrections, but it is plausible to deduce that crit- 
ical exponents are converging to the Ising EA results, 
as the experimental results for the IIB Mn Te suggest 
Further open questions need clarification: a bet- 
ter measure of the critical exponents using the method 
of the precise characterization of the AFM phase, 
and a detailed study of the order of the paramagnetic 
anti-ferromagnetic phase transition at low dilution. 

We acknowledge interesting discussions with D. Brogi- 
oli, B. Coluzzi, A. Geddo-Lchmann, G. Parisi, F. Ricci- 
Tersenghi, J. J. Ruiz-Lorenzo. V.M.M. is a M.E.C. fellow 
and has been partially supported by CICyT(AEN97-1708 
and AEN99-1693). The simulations have been performed 
using the Pentium clusters of the Universita di Cagliari 
(Kalix2) and the Universidad de Zaragoza (RTNN col- 
laboration) . 
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1.277(22) 


1.91(7) 
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TABLE I. Transient critical exponents u{Li,L2) and 
rj{Li,L2) obtained with the quotient method (see text). 
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